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We study various measures of classicality of the states of open quantum systems subject to 
decoherence. Classical states are expected to be stable in spite of decoherence, and are thought 
to leave conspicuous imprints on the environment. Here these expected features of environment- 
induced superselection (einselection) are quantified using four different criteria: predictability sieve 
(which selects states that produce least entropy), purification time (which looks for states that are 
the easiest to find out from the imprint they leave on the environment), efficiency threshold (which 
finds states that can be deduced from measurements on a smallest fraction of the environment), and 
purity loss time (that looks for states for which it takes the longest to lose a set fraction of their initial 
purity). We show that when pointer states - the most predictable states of an open quantum system 
selected by the predictability sieve - are well defined, all four criteria agree that they are indeed the 
most classical states. We illustrate this with two examples: an underdamped harmonic oscillator, for 
which coherent states are unanimously chosen by all criteria, and a free particle undergoing quantum 
Brownian motion, for which most criteria select almost identical Gaussian states (although, in this 
case, predictability sieve does not select well defined pointer states.) 

PACS numbers: 03.65.Yz, 42.50.Lc 



I. INTRODUCTION 

Persistent monitoring of an open quantum system by 
its environment can single out a preferred set of states, 
known as pointer states. Pointer states are the most 
robust in spite of the interaction with the environment. 
That it, they entangle least with the environment, and, 
hence, are least perturbed by it 0,0- This is the essence 
of environment-induced superselection or einselection. 

In the standard treatment of environment-induced de- 
coherence, the search for pointer states is based on the 
unconditional master equation for the reduced density 
matrix of the system [j. [J. 0]. Unconditional dynam- 
ics is obtained after tracing over the environment, i.e., 
by discarding any information about the system present 
in the environment. Predictability sieve, introduced in 
0, H, explored and extended in 0, and reviewed 
in 0, S @' i s a wa y ^0 quantify predictability and, 
hence, classicality of states: For every pure initial state 
I^P) of the system one asseses the loss of predictabil- 
ity caused by the environment by computing entropy 
H\ty) = — Trp|^ (t) logp^) (<), or some other measure of 
predictability (e.g., purity Tr(p|^) 2 ) from the reduced 
density matrix P|^)(i) (which starts as p\\b)(0) = \^f)(^\). 
Entropy is a function of time and a functional of the ini- 
tial state 1^). Pointer states are obtained by minimizing 
H\q,j over l^) and demanding that the answer be robust 
when varying t within a reasonable range. 

For example, in an underdamped harmonic oscillator 
coherent states are the pointer states . This means that 
for any time longer than a fraction of the period of the 
oscillator, starting the evolution from a coherent state 
(or some slightly squeezed state that is very close to a 
coherent state) maximizes predictability. For very short 



times (i.e., times short compared to the oscillator pe- 
riod) states squeezed in position tend to do better when 
the oscillator is coupled to the environment through posi- 
tion. These subtleties disappear when the rotating wave 
approximation is adopted, and we shall do so here. But, 
generally, for such short times (or, alternatively, for a free 
particle) equally unambiguous pointer states do not exist: 
the most predictable states tend to depend significantly 
on the time t for which the prediction is needed. 

Recently, it has been recognized that one finds out 
states of systems indirectly by measuring the environ- 
ment [j| ■ This recognition elevates the role of the envi- 
ronment to a communication channel and to a witness 
of the state of the system 0] . The dynamics of this pro- 
cess can be analyzed using unravellings of the uncondi- 
tional master equation for the density matrix of the sys- 
tem [HEH G3, El El El This is not the only possible 
approach: information theory can be employed to anal- 
yse the role of the environment as a witness perhaps even 
more rigorously [UE,E,E]. 

The program of quantum unravellings was initiated 
some time ago and developed primarily in the context 
of quantum optics [l9j . More recently, it has been ap- 
plied to study the problem of decoherence and the quan- 
tum/classical transition [H O, E E Q QS HJ . For a 
fixed system-environment interaction, different quantum 
unravellings correspond to different measurements on the 
environment. 

Our work [T3. suggested that pointer states should 
be the easiest to find out from the environment: The 
same measurement scheme whose records correlate with 
the pointer states is the one that has the shortest pu- 
rification time of any initial mixed state. This, as well 
as related results [l7|, lead one to conjecture that differ- 
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ent classicality criteria (such as predictability sieve and 
purification time) may lead to the same pointer states. 
This conclusion is also of obvious interest in applications 
(e.g., quantum control) where one would like to tailor 
the measurement strategy to optimize some criteria of 
performance. 

In this paper we will analyze further different clas- 
sicality criteria, some of them considered recently in 
Refs . fl5j . and we will show that when pointer states 
are well defined, all these criteria are indeed optimized 
by a single strategy for monitoring the environment (that 
always selects the same pointer states of the system). 
When pointer states are imperfect, different criteria are 
optimized by different measurement schemes. However, 
at least in the examples we have encountered, all these 
schemes still single out almost identical, most classical 
states. 

We consider two models. We first study an under- 
damped harmonic oscillator coupled to a zero tempera- 
ture and to a finite temperature harmonic oscillator bath. 
As noted earlier, in this case coherent states are pointer 
states. We shall show that a single unravelling - het- 
erodyne detection - optimizes all four classicality cri- 
teria that we consider. We also study the model of a 
free particle undergoing quantum Brownian motion at fi- 
nite temperature. A free particle does not have equally 
well defined pointer states: they are Gaussians that are 
approximately position eigenstates at high temperature, 
and less squeezed Gaussians at low temperatures, but 
the exact "shape" of the Gaussians depends on t. We 
show that, although in this case different classicality cri- 
teria are optimized by different measurement strategies, 
the states that are singled out are nearly identical in all 
cases. 



II. AN UNDERDAMPED HARMONIC 
OSCILLATOR 



describes damping and decohcrcncc due to spontaneous 
emission of phonons, and 

H4a]p = ape^ + paV^-pTrlape-^ + paV 1 *], (3) 

is a nonlinear "discovery" term that feeds back into the 
master equation information about the state of the sys- 
tem gained from measuring the environment. We use 
units where the spontaneous emission time is 1. Here 
V = Vx+Vy is the measurement efficiency, 4> is the phase of 
the local homodyne oscillator, dW x and dW y are uncor- 

0, 



related Gaussian Wiener increments, dW T 



dW„ 



dWl = dWg = dt, and dW x dW y = 0. The choice t] x = r] 
and rj v = corresponds to homodyne unravelling, while 
Vx = Vy — i 7 ? to heterodyne detection 0,H2- As the 
harmonic oscillator Hamiltonian H is invariant under ro- 
tations in the x — p plane from now on we set (f> = 
without loss of generality. After fixing (f> the measure- 
ment scheme is fully determined by the pair of numbers 

It is useful to transform Eq.Q to an equation for the 
Wigncr probability distribution W(x,p): 

dW(x,p) = dtDW + ^dW x [HoW-WTr(H W)} + 
^dWylH^W-WTriH-W)], (4) 

where 

1 , „ „ s 1 



VW=[1 + -(xd x +pd p ) + -{di + flg) ) W(x,p)dt, 



(5) 



and 



H^W = V2 (x cos (f) + p sin (ft+ 

^(coscf)d x +sin<f)d p )^W(x,p). (6) 

Since the unconditional (77 = 0) steady-state has a Gaus- 
sian Wigner distribution, and the conditional evolution 
preserves Gaussianity, we propose a solution of the form 



A. Zero temperature bath 

We consider a one dimensional harmonic oscillator 
with frequency w, position x and momentum p in a zero 
temperature environment of harmonic oscillators linearly 
coupled to the system oscillator. We consider continuous 
Markovian unravellings with homodyne and heterodyne 
detection schemes for the phonon environment. In the 
interaction picture and after the rotating wave approx- 
imation, which assumes that to is much larger than the 
spontaneous emission rate, the stochastic master equa- 
tion (SME) describing the conditional evolution of the 
density matrix of the system is 

dp = V[a]pdt+^dW x H^[a]p+^dW y H^[a]p, (1) 

where the Lindblad term 



T>[a]p = apa) — —a^ap 



1 



pa' a, 



(2) 



W(x,p) = e 



-a(t)[x-x (t)] 2 -/3(t)[p-p (t)] 2 +S(t) 



(7) 



It is also worth noting that there exists a general argu- 
ment that the most predictable states ought to be Gaus- 
sians 0, 0. The inverse of the variances in position, 
a(t), and in momentum, (3(t), evolve deterministically 
even under the conditional stochastic evolution: 



a(t) = 

m = 



e t [a {l - rj x ) + r] x ] - (1 - a )r] x 
e*[a (l - Vx) + Vx) + (1 - «o)(l - Vx) 
e t [/3o(l-%)+ % ]-(l-/3 )% 



. (8) 



e t [/3o(l-%)+^] + (l-/3o)(l-r / j / )' 

where «o and /3o are the initial conditions. In particular, 
the unconditional evolution with r\ x = r\ v = 0, 

f 



a(t) 



1 - e-*(l - l/a ) 
1 



(9) 
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FIG. 1: Purity as a function of time under unconditional 
evolution at T = starting from an initial pure Gaussian state 
parametrized by a squeezing parameter £ = log k. Purity 
is not lost for £ = 0, i.e., the coherent state with (a,/3) — 
(1, 1). Coherent states are perfect pointer states at T = 0, and 
optimize the predictability sieve criterion. By contrast, purity 
is initially lost on a decoherence time scale for other states 
with £ 7^ 0. Time is measured in units of the spontaneous 
emission time. 



relaxes the state to vacuum on a timcscalc t ~ 
max [ln(l/ao)j m (l//3o)]- F° r an initial thermal state 
with ao = fto << 1 this relaxation time is much larger 
than 1. 

The first criterion we consider is predictability sieve. It 
is well known that for an underdamped harmonic oscilla- 
tor at zero temperature, preferred states - as selected by 
the predictability sieve - are coherent states 0, : these 
are the initial pure states that minimize the rate of purity 
loss. In fact, when T = they are perfect pointer states 
in the sense that an initial coherent state (ao = /?o = 1) 
remains a pure coherent state with a(t) = 0(t) = 1 , se e 
Eqs.JSJ), and the purity loss is identically zero |23 . l24ll25l| . 
In Fig. Q we plot purity P(t) = y?a(t)[3(t) for zero 
temperature under unconditional evolution starting from 
initial pure states parametrized by one parameter k as 
a(0) = k, /3(0) = 1/k. As already mentioned, the pre- 
dictability sieve criterion is optimized for initial coherent 
states (k = 1). Now we will show that these states are 
also the most classical states according to three differ- 
ent criteria: purification time, efficiency threshold, and 
purity loss time. We will also show that a single mea- 
surement strategy, heterodyne detection, maximizes the 
robustness of all those criteria. This is expected, since 
heterodyne detection measures the x and p quadratures 
of the state on an equal footing, and a coherent state has 
equal variances of position and momentum. 

The second criterion we consider is purification time. 
We start from a thermal state with ao = /3q -C 1, and 
study how fast the state is purified by the measurement. 
Purity P(t) = -J a(t)f3(t) can be shown to evolve deter- 
ministically, so there is no need to do conditional stochas- 




FIG. 2: Purification time as a function of the measurement 
scheme parametrized by s £ [0, 1]: r/ x — 77 sin 2 (s-|) and r) y = 
7/ cos 2 (sf). The initial state is a fat Gaussian state ao = 
[3o = 10 -8 . The temperature of the bath is zero. We note 
that while there is a clear optimum in Fig. (1), the very 
best state s = 0.5 is not all that much better than s ~ 0.2 
or 0.8. Thus, there is a large "quantum halo" [2r3 | of states 
that maintain purity well (and in this sense are reasonably 
classical) . 



tic simulations to study purification time. The initial rate 
of purity gain can be easily evaluated 



dP 
~dt 



VVxVy- 



(10) 



t=o 



It is maximized for heterodyne detection r\ x = rj y = i]/2. 
To go beyond this simple calculation of the initial pu- 
rification rate, we calculate purification time, defined as 
the time needed by the conditional evolution to increase 
purity to P = 0.5, which is half way between the initial 
thermal value P = 0, and the final pure state P = 1 [Ifij . 
In Fig.(l) we plot purification time for different values of 
the global efficiency 77 as a function of the measurement 
scheme, parametrized by s £ [0,1]: r\ x = Tycos 2 (sf) 
and r\ y = 77 sin 2 As expected, purification time is 

minimal for heterodyne detection (s = 0.5), which corre- 
sponds to a coherent state (i.e., a pointer state). 

The third criterion is efficiency threshold: the mini- 
mal efficiency 7/ t h r needed to reach a threshold value of 
purity < Pthr < 1 starting from a mixed high tempera- 
ture thermal state. This criterion cannot be directly used 
here to discriminate between different measurements be- 
cause in a long time the state always relaxes to the pure 
vacuum state no matter what the efficiencies 77^ and 77^ 
are, even in the unconditional evolution (28|. When ap- 
plied literally, the criterion always gives the trivial result 
Vthr = 0. This is why we study a modified criterion: 
we measure the minimal 77ti ir needed to reach a P t hr af- 
ter a time tthr- In Fig- EH we show ?7thr needed to reach 
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FIG. 3: The minimal efficiency threshold r/thr needed to reach 
the threshold purity P t hr = 0.5 after the time tthr = 2.5 as 
a function of the measurement scheme parametrized by s £ 
[0, 1]: rj x = ?7th r sin 2 (sf ) and % = i] t hv cos 2 (s%). The initial 
state and temperature are the same as in Fig. I. Again, 
s = 0.5 is surrounded by a sizable "halo" of (nearly) pointer 
states I2I. 



-Pthr = 0.5 after i t hr = 2.5 for different measurement 
schemes parametrized by s e [0, 1]: r\ x = ?ycos 2 (s§) and 
T] y = 77 sin 2 (s^). As expected, the r^hr is minimal for 
heterodyne detection (s = 0.5). 

The fourth criterion is the purity loss time (proposed 
in and called "mixing time" in that reference). The 
conditional trajectory has to evolve long enough for the 
state to become pure. Then one stops measuring the 
environment to see how fast the purity of the state de- 
creases as a function of time, for different initial measure- 
ment schemes. All this has to be done before the state 
relaxes to vacuum. As already mentioned, for a high 
temperature initial state with ao = /3q << 1 the uncon- 
ditional evolution relaxes the state to vacuum after time 
ln(l/ao) ^ 1- In the limit of ao — > this relaxation 
time tends to infinity. In the same limit the conditional 
evolution becomes 



a(t) 



e'-l 



e« - 1 

'ly 



(11) 
(12) 



and the purity P(t) = -J a(t)(3(t) approaches 1 after a 
purification time max [hi(l/ri x ), ln(l/r] a )]. When both 
1 > VxiVy ^ a o, then the purification time is much 
shorter than the relaxation time (compare also Fig|2J), 
and the conditional state collapses to a coherent state 
much faster than the relaxation time in the unconditional 
evolution. To measure the purity loss time once the pure 
coherent state is achieved, the evolution is switched to 



the unconditional evolution with r\ x = r\ y = 0, and the 
purity decay is observed. However, since all measure- 
ment strategies collapse the state to a coherent state, the 
subsequent unconditional evolution does not destroy the 
purity of this state: the amplitude of the pure coherent 
state smoothly relaxes to the vacuum. This infinite pu- 
rity loss time is the same for all measurement schemes. 
Coherent states, in the case of T = 0, are perfect pointer 
states: once the conditional state is collapsed to a coher- 
ent state it cannot lose any purity. 



B. Finite temperature bath 

The SME for the harmonic oscillator coupled to a finite 
temperature bath at temperature T is 

dp = (n + l)D[a]pdt + nD[a)]pdt + 



Vl + 2n 



dW x {H4,[(n + l)a]p - H^nc 
dW v {H^[{n + l)a]p 



T~ta 



\ [na f ]p} , 



where n = (e fiw °/ feBT — is the Bose distribution. This 
equation is valid in the interaction picture and after the 
rotating wave approximation is performed. We can write 
the analogous equation for the Wigner distribution. We 
take — 0, without loss of generality, and get 

dW(x,p) = dtDW + ^dW x [H x W-WTr(H x W)} + 
^dWy[H y W-WTr(H y W)], (13) 

where now 



and 



A + xd x +pd z + l_ + 2n 2 2 



X 1 p ■ 



W(x,p)dt, 



H X W = V2(x+^(l + 2n)d x jW(x,p), 
H y W = V2\p+^(l + 2n)d p jW(x,p). 

Assuming a Gaussian Ansatz for the Wigner function we 
get the deterministic equations for a(t) and (3{t). They 
are the same equations as for zero temperature, with the 
replacement a(t) — > (1 + 2n)a(t). Hence, the solution 
can be straightforwardly obtained from Eq.(JH|. 

When T > 0, the stationary state of the evolution is 
a ss = f3 ss = 1/(1 + 2n), i.e., a Gaussian with variances 
in position and momentum larger than 1. The uncondi- 
tional evolution relaxes to this stationary state on a time 
scale t ~ max[ln(l/(l + 2n)ao),ln(l/(l + 2n)/3 )]. An 
initial pure coherent state ao = /3q = 1 does not remain 
pure; it loses purity and reaches the stationary state with 
final purity P(t -> 00) = 1/(1 + 2n) < 1. Predictabil- 
ity sieve still selects coherent states as pointer states, 
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FIG. 4: Purity as a function of time under unconditional 
evolution at finite temperature starting from an initial pure 
Gaussian state parametrized by the squeezing parameter £ = 
logs. The temperature of the bath is k-sT/fiuio — 10 6 , where 
ujo is the frequency of the system oscillator. Again, coherent 
states (£ = 0) are selected by predictability sieve, as they lose 
purity most slowly. For large times, all initial states lead to 
the same stationary state with purity P(oo) — 1/(1 + 2n), 
where n « 10 6 for this temperature. Time is measured in 
units of the spontaneous emission time. 



even at finite temperature (see Fig. Q). Also, heterodyne 
detection measurement maximizes the robustness of all 
the criteria (purification time, efficiency threshold, and 
purity loss time). For an initial high temperature state 
such that (1 + 2n)ao = (1 + 2n)/3 <C 1, the uncondi- 
tional evolution relaxes towards the stationary state on 
a very large time scale with P(oo) = 1/(1 + 2n), while 
the conditional evolution reaches the same asymptotic 
value of purity on a much shorter time scale, of the order 
of max[ln(l/?7 x ), ]n(l/r]y)]. At finite temperature, the pu- 
rity loss time is not infinity as for T = 0: all measurement 
schemes lead to the same, non pure state, then the mea- 
surement is switched off, and then purity is reduced again 
to the same final value, irrespective of the measurement 
scheme. This criterion cannot distinguish measurement 
schemes. 

In summary, all classicality robustness criteria select 
the same states (coherent states) for the model of an 
undcrdamped harmonic oscillator, for all temperatures 
of the bath. 



III. A FREE PARTICLE UNDERGOING 
QUANTUM BROWNIAN MOTION AT HIGH 
TEMPERATURE 



The high temperature conditional homodync and het- 
erodyne master equation for a free particle undergoing 



dp 



qp + pq 



dt 



[y/rj dW (cp-pTr(pc)) 



>rV[c]pdt' 

h.c] . 



(14) 



Here r\ is the total efficiency of the measurement, and the 
annihilation operator c is 



(15) 



We use rescaled units such that the damping rate, the 
particle mass, Boltzmann constant, and h are all unity. 
The complex, Gaussian Wigner increment has correlators 
dWdW* = dt and dW 2 ' 




dt 



ith -1 < r < 1. It 



can be written in terms of two real, uncorrelated, Gaus- 



sian noises dW x and dW y as 



dW 



1 + r 



dW T 



,-i<p/2 



1 - r 



dW y , (16) 



where dW£ = dWy 2 = dt, and dW x dW y = 0. In this way, 
the stochastic part of the master equation Ea. H14|) can 
be expressed as 



dp' 



stoch 



Rf x dW x H^cjp + ^/%dW y H 0+f [c]p. (17) 



Here we have defined <f> = —ip/2, r\ x = rj(l + r)/2, and 
Vy = ^(1 — r )/2- Heterodyne detection corresponds to 
r = 0, while r = 1 corresponds to homodyne measure- 
ment of the linear combination ce 1 * + c^e -4 ^ of q and p. 
Hence, = corresponds to a measurement of q, and 
cf> = 7i"/2 to a measurement of p. The unconditional ver- 
sion of the master equation fTH) was derived in Ref.j23 as 
Markovian approximation to the exact quantum Brown- 
ian motion master equation at high temperature. The 
equation is valid only when T>1. 

As discussed in the Introduction, the free particle un- 
dergoing quantum Brownian motion does not have well 
defined pointer states for any temperature of the envi- 
ronment. For example, at high T the Lindblad term in 
the master equation i|14|) formally dominates over the 
free Hamiltonian part, and predictability sieve would se- 
lect eigenstates of the annihilation operator c, that corre- 
spond to Gaussian states squeezed in position. However, 
this answer is not robust under variations of the time 
t during which the criterion is applied: states squeezed 
in position have a large dispersion in momentum, and 
the free dynamical evolution will drive the state away 
from the squeezed eigenstates of c. We will show below 
that, because of these facts, different classicality criteria 
are maximized for different unravellings. However, even 
under these circumstances the selected states are very 
similar Gaussians. 

Just as in the model of the underdamped harmonic 
oscillator, the conditional dynamics Ea. (|14f> preserves 
Gaussianity. We assume a Gaussian ansatz for the 
Wigner function 

W(x,p) = cxp{-a{t)[x-x Q {t)] 2 ~P(t)[p-p Q (t)} 2 - 
2 1 (t)[x ~ x (t)][P - Po(t)] + S(t)} . (18) 
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The coefficients a(t), f3(t), and 7(f) evolve dcterministi- 
cally according to the following equations 



da 
~dt 
d£ 
~di 

~dt 
where 



4T 

-ATP 2 
or) 



4T 7 2 + ATr, x A\ + AT Vy A 2 2 , (19) 
2 



rjy_ 
AT' 



D 



2 > 



= - ^T[3 1 + 7 - a + r) x A x B x - 7 ly A 2 B 2 , 



7 sine 



1 



AT) 



COS ( 



a \ 

1 sin ( 

AT) 



A 2 = 7Cos< 

Bi = 7COS0- (1 - 4T/3)sin0, 
B 2 = 7sin0+ (1 - 4T/3)cos0. 

To apply the predictability sieve criterion we solve the 
unconditional version of these equations (rj x = rj y = 0) 
starting from different initial pure Gaussian states. The 
initial pure state is parametrized by two parameters 
(A,C) as a(0) = 4T A, 7(0) = C and /3(0) = [1 + 
7 2 (0)]/a(0). In this parametrization the squeezed eigen- 
state of c is (A = 1, C = 0). In figure El we show purity 
P(t) = y/a(t)/3(t) — 7 2 (£) after unconditional evolution 
for a very short time t = 0.1/4T = 2.5 10~ 8 at temper- 
ature T = 10 6 as a function of initial pure state (A, C). 
Purity is maximal when (A, C) are approximately equal 
to (1, 0), i.e., purity is maximized for the eigenstate of the 
squeezed operator c with (a, (3, 7) = (4T, 1/4T, 0). As ex- 
pected, at early times the evolution is dominated by the 
Lindblad term and it is most predictable for the eigen- 
states of c. In the high-T limit these eigenstates tend 
to be position eigenstates. However, the range of values 
of purity in figurcElis only (0.99999999, 1), so that after 
this short time the evolution is still very predictable no 
matter what the initial state is - the most robust states 
are not well distinguished from the rest. 

Rather than the early times t ~ 1/4T we found the 
times t ~ 1/V4T more relevant from the point of view 
of predictability sieve, since then the unconditional evo- 
lution starting from the most robust initial states loses 
approximately one half of its initial purity, while at the 
same time the evolution starting from bad states is es- 
sentially not predictable at all. At t ~ \j\f4T the most 
robust states are well distinguished from "the Hilbcrt 
space chaff' by the predictability sieve. In figure we 
show purity P(t) after unconditional evolution at tem- 
perature T = 10 6 at a time t = 2/VAT = 1(T 3 . Pu- 
rity is maximal when (A, C) = (1.75,1.75) or (a,/3, 7) = 
(1.75^4^,^,1.75). In the high T limit these states 
are very different from the eigenstates of c although 
they also tend to approximate position eigenstates for 
large values of the temperature. In a strict sense pre- 
dictability sieve does not select any well defined and t- 
independent pointer states. Nevertheless, non-local su- 
perpositions are still quickly destroyed, and lose purity 
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FIG. 5: Purity under unconditional evolution at T = 10 6 
after a short time t = 0.1/4T = 2.5 10 -8 starting from an ini- 
tial pure Gaussian state parametrized by two numbers (^4, C). 
The parameter A is related to the squeezing of the prob- 
ability distribution of the initial state, and C is related to 
its tilt in phase space. Purity is maximized for an initial 
state that is approximately parametrized by (A,C) = (1,0), 
that is, the eigenstate of the squeezed annihilation operator 
c with (q,/3, 7) = (4T, j=,0), which is the most predictable 
state after this short time of unconditional evolution. How- 
ever, the vertical range of the plot is only (0.99999999, 1) - at 
this early time the unconditional evolution is still very pre- 
dictable no matter what its initial pure state is. The most 
predictable eigenstates of c are not well distinguished from 
the rest. Rescaled units are used. 



much faster than the stable and more or less localized 
Gaussians we have described. Therefore, for many prac- 
tical purposes the states which are selected at measure- 
ment times t ~ l/v4T may play the role of the robust 
classical states. Parameters of these states scale with T 
as a ~ V4T, /3 ~ yjf > an< ^ 7—1- In the following we 
will see that this class of states is also preferred by the 
other classicality criteria. 

Now we turn to the second classicality criterion, pu- 
rification time. In order to quantify it we prepare the 
system in its unconditional stationary thermal state, 
a a = 7o = 0,/3o = 1/2T. Then we compute how fast 
the conditional state is gaining purity P(t) depending on 
the measurement scheme on the environment, defined by 
r and <p. Here we assume full efficiency (77 = 1). The 
initial rate of gain of purity can be easily calculated an- 
alytically at t = 0: 



1 + r cos 2(f> 



and it is found not to depend on T. This initial pu- 
rification rate is fastest when r = 1 and cj> = 0, tt, or 
when r = — 1 and <f> = ir/2, 3tt/2. All these four solutions 
correspond, in fact, to the same x— homodyne. The slow- 
est initial purification is obtained for y— homodyne with 
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purity 




FIG. 6: Purity after unconditional evolution at T = 10 for 
a time t = 2/\/4T = 10~ 3 when the unconditional evolution 
starting from the most robust initial state loses approximately 
one half of its initial purity. Purity is maximal when (A, C) = 
(1.75, 1.75), or (a,0,j) = (1.75-v/lT, ^§=,1.75). At this time 
the most predictable initial states are well distinguished from 
(are much more pure than) other states, which already have 
negligible purity. Rescaled units are used. 
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homodyne angle 



FIG. 7: Purification time as a function of the homodyne 
angle cj> ( r = 1) f° r an initial thermal state with temperature 
T. Purification time is the shortest for x— homodyne, i.e., 
when (f> — 0, n. By contrast, y— homodyne (<f> = 7r/2) does 
not purify the state at all. Rescaled units are used. 



r = 1, <j) = 7r/2. To go beyond the initial purification we 
measure the purification time, defined as the time needed 
by the fully efficient (rj = 1) conditional evolution to in- 
crease purity to P — 0.5, which is half way between the 
initial thermal state value (P as 0) and the final pure 
state value (P = 1) 15|. Purification times for different 
homodyne angles (f> a t different temperatures are shown 
in figure [7| As expected, purification time is the short- 
est for x— homodyne, i.e., when r = 1 and <f> = 0, 7r. In 
contrast, y— homodyne does not purify the state at all. 

Now we consider the efficiency threshold, defined as the 
minimal efficiency 7jthr required for the conditional evo- 
lution to asymptotically reach a threshold purity Pthr- 
Here we assume P t \ u - = 0.5 and an initial thermal state. 
Efficiency thresholds for different homodyne angles <fi 
are collected in Figure OH The threshold is the low- 
est for x— homodyne ((f) = 0, n), and the highest for 
y— homodyne when full efficiency is required to find out 
the state after infinitely long time. 

Finally, in order to measure the purity loss time (called 
"mixing time" in ^H)' we evolve the conditional state of 
the system for a sufficiently long time so that it reaches 
a conditional stationary state, that depends on the mea- 
surement scheme. The stationary conditional state can 
be obtained by setting the left hand sides of Eqs. (|19|) to 
zero and solving the equations for a ss , /3 SS , 7 SS . In the high 
T limit, the stationary solutions can be found in the form 
a ss = A SS V4T,P SS = B ss /V4T, 7ss = C ss . After keeping 
in each equation only the leading order terms in the high 
T limit we obtain simplified equations for A ss , B SS ,C SS , 



2 0.8 
o 

J 0.6 

>, 

o 

g0.4 



a3 0.2 



vT= 1000000. 



T= 1 
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FIG. 8: Efficiency threshold as a function of the homodyne 
angle <j> for an initial thermal state with temperature T. The 
threshold is the lowest for x— homodyne {cj> = 0, n) and the 
highest for y— homodyne {cj> — tt/2). Rescaled units are used. 



and their solution gives, 

-^ss Css 



^4ss — 
Bss = 



(1 + r cos 2(f)) — r sin 2(f), 



4C SI 



1 + r cos 2(f> 



1/2 



r sin 2(f) + + 2r cos 20 + r 2 

^ss — 7— ■ ~ L • [2U) 

1 + r cos 2(p 

Once the stationary state is achieved the evolution is 
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FIG. 9: Purity loss time needed by unconditional evolution to 
reduce purity from the initial value P = 1 to P = 0.5, which 
is half way to the stationary value P = in the thermal state. 
The homodyne angle that maximizes the purity loss time is 
cj> — 1.35 ^. On the other hand, the purity loss time is zero 
for y— homodyne with (j> = 7r/2. Rescaled units are used. 



switched to unconditional evolution with r\ = and 
we calculate the purity loss rate. This depends on the 
measurement scheme through the stationary conditional 
state, which is the initial state for the unconditional evo- 
lution. The initial purity loss rate can be easily calculated 
as 



dP 
~~dt 



= Vtb s: 



(21) 



to leading order in high T . This initial purity loss rate is 
minimized for a;— homodyne with r = 1 and homodyne 
angle <p = 

So far we have analyzed only the initial purity loss 
time. Now we consider the purity loss time, defined as 
the time needed by the unconditional evolution to reduce 
purity to P = 0.5, which is half way between the initial 
pure state and the stationary thermal state with PsO. 
In figure [5] we show the purity loss time as a function 
of homodyne angle (f> after homodyne measurement with 
r = 1. The homodyne angle that maximizes the purity 
loss time is <fi — 1.357r/2. 

We applied the three different criteria and we found 
that in two cases - purification time and efficiency 
threshold - the most robust measurement strategy is 
x— homodyne ((/> = 0) over the whole range of tempera- 
tures, but in one case - purity loss time - the most robust 
measurement is homodyne with <j) = 1.357r/2. In each 
case the most robust measurement defines the most clas- 
sical state as the conditional stationary state obtained 
with the optimal measurement scheme. These results may 
leave one under the impression that various criteria lead 
to different candidates for classical states. Indeed, this 
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100 



FIG. 10: Wigner functions of the stationary conditional states 
selected at high temperatures T = 10 6 by different robustness 
criteria: (a) purification time and efficiency threshold (homo- 
dyne angle <j> = 0) , and (b) purity loss time (homodyne angle 
cj> = 1.357r/2). The two selected states are very similar (scalar 
product equal to 0.97), and Gaussians squeezed to "approxi- 
mate" position eigenstates. Rescaled units are used. 



was the conclusion of [l5j , where similar results concern- 
ing optimal measurement strategics were obtained for the 
model of this Section. The key point is, however, that 
although different classicality criteria are maximized by 
different measurement strategies, in all cases the states 
that are singled out are almost identical. 

In the two following tables we list parameters of the 
conditional stationary states for <fi = and <f> = 1.357r/2 
for several values of the temperature. The states can 
be described by three independent quantities: the dis- 
persions in position Aa; and in momentum Ap, and the 
covariance C xp = (xp + px) /2 — (x) (p). We find that for 
all listed temperatures the stationary conditional states 
for <j) — and <j) = 1.357r/2 are very similar (see figures 6 
and 7). One can easily imagine they are well within each 
others "quantum halo" [2(| . These states evolve from al- 
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FIG. 11: Wigner functions of the stationary conditional states 
selected at low temperatures T — 1 by different robustness 
criteria: (c) purification time and efficiency threshold (homo- 
dyne angle <j> = 0), and (d) purity loss time (homodyne angle 
<j> — 1.357r/2). The two selected states are very similar (scalar 
product equal to 0.99), and approximately equal to coherent 
states with a = /3 = l,7 = (scalar product with a coherent 
state is equal to 0.99 and 0.98 respectively). Rescaled units 
are used. 



most coherent at T = 1 to almost position eigenstates at 
high T. Thus, in contrast to ^5|, we conclude that all 
of the above criteria lead to essentially the same "most 
classical" states. 



TABLE I: Stationary conditional state for homodyne angle 
4> — at different temperatures T. Rescaled units are used. 



T 


a S s 


f3 ss 


7ss 


Ax 


Ap 


Cxp 


10 b 


2826 


0.0007 


-0.999 


0.018 


37.968 


0.509 


10 4 


280 


0.0070 


-0.992 


0.059 


11.977 


0.508 


10 2 


26.4 


0.0707 


-0.931 


0.188 


3.633 


0.466 


1 


1.53 


0.8002 


-0.480 


0.632 


0.877 


0.241 



TABLE II: Stationary conditional state for homodyne angle 
cj> = 1.357T-/2 at different temperatures T. Rescaled units are 
used. 



T 




A s 




Ax 


Ap 


C X p 


10 b 


1500 


0.0007 


-0.281 


0.018 


27.792 


0.144 


10 4 


149 


0.0072 


-0.280 


0.060 


8.656 


0.141 


10 2 


14 


0.0741 


-0.270 


0.192 


2.694 


0.140 


1 


1.05 


0.9561 


-0.112 


0.694 


0.728 


0.056 



IV. CONCLUSIONS 

In the example of the under-damped harmonic oscil- 
lator, where pointer states defined by the predictability 
sieve criterion are unambiguous, we have shown that dif- 
ferent classicality criteria single out the same states, i.e., 
coherent states that coincide with the pointer states se- 
lected by the predictability sieve @, 3 . We considered 
three of these "alternative classicality criteria" : purifica- 
tion time, which looks for states of the system that are 
the easiest to find out from the imprint they leave on 
the environment; efficiency threshold, which finds states 
of the system that can be deduced from measurements 
on a smallest fraction of the environment; and purity 
loss time, that looks for states of the system for which 
it takes the longest to lose a set fraction of their initial 
purity. Our findings indicate that quantum unravelings 
are in effect classical for pointer states. 

When pointer states are not well defined, as in the 
model of a free particle undergoing quantum Brownian 
motion, or when pointer states do not exist at all, dif- 
ferent criteria may select different states. However, it 
appears that these candidate classical states are very sim- 
ilar, defining approximate pointer states, as in the model 
of the free particle undergoing quantum Brownian mo- 
tion. This is an interesting example because different 
classicality criteria select measurement schemes that ap- 
pear quite different, but these different schemes turn out 
to prepare very similar Gaussian states. 

Following our work 0] > the classical robustness of dif- 
ferent quantum unravcllings was studied in |l5j . where 
some of the classicality measures used in this paper - such 
as efficiency threshold and purity loss time (called "mix- 
ing time" in [l^) - were considered. It was claimed there 
that for a fixed environmental interaction the level of ro- 
bustness depends on the measurement strategy, and that 
no single strategy is maximally robust in all ways. This 
conclusion was drawn from two models: resonance fluo- 
rescence of a two-level atom (a model for which pointer 
states do not exist), and a free particle in quantum Brow- 
nian motion (one of the models used in this paper, for 
which one might perhaps argue there are approximate, 
time-dependent, and imperfect pointer states). As we 
discussed above, the conclusion of Ref.jl^ does not ap- 
ply to the case when einselection works well, i.e., when 
pointer states arc well defined (e.g., the under-damped 
harmonic oscillator). Moreover, the conclusion of [l5j 
is misleading even for the example of the free particle 
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undergoing quantum Brownian motion they have inves- 
tigated: Pointer states are not well defined there, and 
different robustness criteria are optimized by different 
measurement strategies. Yet, the states resulting from 
all of the above strategies are essentially identical! Thus, 
even when einselection does not pick out unique pointer 
states, different criteria still agree on what appears to be 
most classical. 

Our conclusions may assist in the choice of the op- 
timal measurement strategies, especially in applications 
that involve quantum control. The fact that when good 
or even approximate pointer states exist, they tend to be 



prepared by all the seemingly different schemes optimized 
for several reasonable criteria, attest to the practical im- 
plications of their classicality. 
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